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Abstract 

This paper focuses on multirate IEEE 802.11 Wireless LAN employing the mandatory Distributed Coordination Function 
(DCF) option. Its aim is threefold. Upon starting from the multi-dimensional Markovian state transition model proposed by 
Malone et.al. for characterizing the behavior of the IEEE 802.11 protocol at the Medium Access Control layer, it presents an 
extension accounting for packet transmission failures due to channel errors. Second, it establishes the conditions under which a 

( s) 

network constituted by A'^ stations, each station transmitting with its own bit rate, , and packet rate. As, can be assumed 
loaded. Finally, it proposes a modified Proportional Fairness (PF) criterion, suitable for mitigating the rate anomaly problem of 
multirate loaded IEEE 802.11 Wireless LANs, employing the mandatory DCF option. Compared to the widely adopted assumption 
of saturated network, the proposed fairness criterion can be applied to general loaded networks. 

The throughput allocation resulting from the proposed algorithm is able to greatly increase the aggregate throughput of the 
DCF, while ensuring fairness levels among the stations of the same order as the ones guaranteed by the classical PF criterion. 

Simulation results are presented for some sample scenarios, confirming the effectiveness of the proposed criterion for optimized 
throughput allocation. 

Index Terms 

DCF, Distributed Coordination Function, fairness, IEEE 802.11, MAC, multirate, non-saturated, proportional fairness, rate 
adaptation, saturation, throughput, traffic, unloaded, unsaturated. 

I. Introduction 

Consider the IEEE802.1 1 Medium Access Control (MAC) layer [1] employing the DCF based on the Carrier Sense Multiple 
Access Collision Avoidance CSMA/CA access method. The scenario envisaged in this work considers N contending stations; 
each station generates data packets with constant rate As by employing a bit rate, , which depends on the channel quality 
experienced. In this scenario, it is known that the DCF is affected by the so-called performance anomaly problem [2]: in 
multirate networks the aggregate throughput is strongly influenced by that of the slowest contending station. 

After the landmark work by Bianchi [3], who provided an analysis of the saturation throughput of the basic 802.11 protocol 
assuming a two dimensional Markov model at the MAC layer, many papers have addressed almost any facet of the behaviour 
of DCF in a variety of traffic loads and channel transmission conditions. 

Contributions proposed in the literature so far can be classified in two main classes, namely DCF Modelling and DCF 
Throughput and Fairness Optimization. 
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DCF modelling. This is the topic that received the most attention in the Uterature since the work by Bianchi [3]. Papers [4]-[6] 
model the influence of real channel conditions on the throughput of the DCF operating in saturated traffic conditions, wlule 
[7]-[9] thoroughly analyze the influence of capture on the throughput of wireless transmission systems. Paper [10] investigates 
the saturation throughput of IEEE 802.11 in presence of non ideal transmission channel and capture effects. The behavior of 
the DCF of IEEE 802.11 WLANs in unsaturated traffic conditions has been analyzed in [11]-[18]. In [19], the authors look 
at the impact of channel induced errors and of the received Signal-to-Noise Ratio (SNR) on the achievable throughput in a 
system with rate adaptation, whereby the transmission rate of the temiinal is modified depending on either direct or indirect 
measurements of the Unk quality. 

Multirate modeUng of the DCF has received some attention quite recently [20]-[24] as well. In [20] an analytical framework 
for analyzing the link delay of multirate networks is provided. In [21]-[22], authors provide DCF models for finite load sources 
with multirate capabiUties, while in [23]-[24] a DCF model for networks with multirate stations is provided and the saturation 
throughput is derived. Remedies to performance anomaUes are also discussed. In both previous works, packet errors are only 
due to colUsions among the contending stations. 

DCF throughput and fairness optimization. This is perhaps the issue most closely related to the problem dealt with in this paper. 
The main reason for optimizing the throughput aUocation of the 802.1 1 DCF is the behaviour of the basic DCF in heterogeneous 
conditions, with stations transmitting at multiple rates: the same throughput is reserved to any contending station irrespective of 
its bit rate, with the undesired consequence that lowest bit rate stations occupy the channel for most time with respect to high 
rate stations [25]. Furthermore, the optimization of the aggregate throughput when different stations contend for the channel 
with different bit rates cannot be done without considering an appropriate fairness approach; the reason is that the optimum 
throughput would be achieved when only the highest rate stations access the channel [25]. In order to face this problem, a 
variety of throughput optimization techniques, which account for fairness issues, have been proposed in the Uterature. Paper 
[25] proposes a proportional fairness throughput aUocation criterion for multirate and saturated IEEE 802.11 DCF by focusing 
on the 802. lie standard. In papers [26]-[29] the authors propose novel fairness criteria, which fall within the class of the 
time-based fairness criterion. Time-based fairness guarantees equal time-share of the channel occupancy irrespective to the 
station bit rate. 

Paper [30] investigates the fairness issue in 802.11 multirate networks by analyzing various time-based fairness criteria. 
It demonstrates that with equal time-share of the channel occupancy among multirate stations, the throughput achieved by a 
reference station in a multirate scenario with A'' contending stations is equal to the throughput that the same reference station 
would achieve in a single rate scenario when contending with other N —1 stations with its same rate. Furthermore, the authors 
prove that the proportional fairness criterion corresponds to fair channel time allocation in a multirate scenario. 

The effect of the contention window size on the performance of the DCF have been also investigated in [31]-[33] in a variety 
of different scenarios. Finally, papers [34]-[39] have been devoted to the throughput optimization of the underlined DCF by 
optimizing a number of key parameters of the DCF, such as the minimum contention window size or the packet size. 

A conmion hypothesis employed in the Uterature regards the saturation assumption, which sometimes does not fit quite 
weU to real network traffic conditions. In real networks, traffic is mostly non- saturated, different stations usuaUy operate with 
different loads, i.e., they have different packet rates, while the transmitting bit rate can also differ among the contending 
stations. Channel conditions are far from being ideal and often packet transmission has to be rescheduled until the data is 
correctly received. Due to Rayleigh and shadow fading conditions, a real scenario presents stations transmitting at different bit 



rates, because of multirate adaptation foreseen at the physical layer of WLAN protocols such as IEEE 802.11b. In all these 
situations the common hypothesis, widely employed in the literature, that all the contending stations have the same probability 
of transmitting in a randomly chosen time slot, does not hold anymore. 

The aim of this paper is to investigate the behaviour of the DCF in the most general scenario of a multirate network, when 
all the previous effects act jointly, as well as to present a proportional fairness criterion which accounts for general loading 
conditions as exemplified by the packet rate As of the contending stations. Contrary to the aforementioned works available in 
the literature, we assume that the s-th station generates data packets with its own size, PL^'^\ with its own constant rate As by 

(s) 

employing a bit rate, , which depends on the channel quality experienced, and it employs a minimum contention window 

(s) 

with size Wq . Moreover, each station is in a proper load condition, which is independent from the loading conditions of the 
other contending stations. Notice that these hypotheses make the model proposed in this work quite different from the ones 
available in the literature, where the saturated condition is mostly adopted. One consequence of the proposed analysis is that 
unloaded, heterogeneous networks do not need any throughput allocation among stations. We propose a theoretical framework 
in order to identify whether a tagged station is saturated, given the traffic conditions of the remaining stations. As a starting 
point for the derivations that follow, we consider the bi-dimensional Markov model proposed in [12], and present the necessary 
modifications in order to deal with multirate stations, non ideal transmission channel conditions, and different packet sizes 
among the contending stations. 

The rest of the paper is organized as follows. Section HIl provides the necessary modifications to the Markov model proposed 
in [12], while the employed traffic model is discussed in Section Hl-DI Section Hill proposes an analytical framework able 
to verify whether a network of N contending stations is loaded. The novel proportional fairness criterion is presented in 
Section IIVI while Section [V] presents simulation results of some sample network scenarios. Finally, Section |VI] draws the 
conclusions. 

II. The Network Scenario: Overview of The Markovian Model Characterizing the DCF 

In [12], the authors derived a bi-dimensional Markov model for characterizing the behavior of the DCF in heterogeneous 
networks, where each station has its own traffic, which could be finite and characterized by the parameter A, expressing the 
packet arrival rate. In order to deal with non-saturated conditions, the traffic model is described by an exponentially distributed 
packet inter-arrival process. In this paper we consider a more general network than [12]. Indeed, in the investigated network, 
each station employs a specific bit rate, R^j^K a specific transmission packet rate. As, transmits packets with size PL^'^\ and it 

Cs) 

employs a minimum contention window with size VFJ , which can differ from the one specified in the IEEE 802.11 standard 
[1] (these modifications are at the very basis of the proportional fairness criterion proposed in Section HVb . A finite retry 
limit is considered in order to avoid infinite number of retries when bad channel conditions inhibit the station from successful 
transmission. 

For the sake of greatly simplifying the evaluation of the expected time slots required by the theoretical derivations that follow, 
we consider Nc < N classes of channel occupancy durationj^- First of all, given the payload lengths and the data rates of 
the N stations, the Nc duration-classes are arranged in order of decreasing durations identified by the index d e {1, • • ■ , Nc}, 
whereby d = 1 identifies the slowest class. Notice that in our setup a station is denoted fast if it has a short channel occupancy. 

'This assumption relies on the observation that in actual networks some stations might transmit data frames presenting the same channel occupancy. As an 
instance, a station STAl transmitting a packet of size 128 bytes at 1 Mbps occupies the channel for the same time of a station STA2 transmitting a packet 
of size 256 bytes at 2 Mbps. 



Furthermore, each station is identified by an index s G {1, • ' ' i ™d it belongs to a unique duration-class. In order to 
identify the class of a station s, we define Nc subsets n{d), each of them containing the indexes of the Ld — \n{d)\ stations 
within n(d), with Ld < N,\/d and J2d=i ^d. — N. As an example, n(3) — {1, 5, 8} means that stations 1, 5, and 8 belong to 
the third duration-class identified by d — 3, and Ld — 3. 

A. Bi-dimensional Contention Markov Model 

The modified bi-dimensional Markov model describing the contention process of the s-th statior|^ in the network is shown 
in Fig. [T] 

Let us elaborate. We consider an overall number of r different backoff stages, starting from the zero-th stage. The maximum 
Contention Window (CW) size is Wmax — 2™Wq'^\ with m < r, whereas the notation Wi = min(2™Wo''\ 2*VKo'*-') is used 

(s) 

to identify the i-th contention window size (Wq is the minimum contention window size of the s-th station). Notice that 
after the m-th stage, the contention window size is fixed to Wmax for the remaining (r — m) stages, after which the packet is 
dropped. An additional backoff stage, identified by (P, — ), with the same window size of the zero-th stage, is considered on 
top of the chain in order to account for the post-backoff stage entered by the station after a successful packet transmission, or 
packet drop [1]. Moreover, the state labelled (P, 0) in Fig. [T] is used for emulating unloaded traffic conditions. 

After the post-backoff stage, a station starts a new transmission because a new packet is available in the queue, provided 
that the channel is sensed idle for DIFS seconds. On the other hand, a new zero-th stage backoff is employed if the channel 
is sensed busy. Notice that the post-backoff stage is entered only if the station has no longer packets to transmit after a packet 
transmission; otherwise, a zero-th stage is started. Moreover, if a new packet arrives during a post-backoff stage, the station 
moves into the zero-th stage, as depicted in Fig. [T] Indeed, backoff stages from to r assume that the station's queue contains 
at least a packet waiting for transmission. 

A packet transmission is attempted only in the states labelled (i, 0), Vi = 0, . . . , r, as well as in the state {P, 0) only if there 
is a packet in the queue and the channel is sensed idle for DIFS seconds. In case of collision, or due to the fact that transmission 
is unsuccessful because of channel errors, the backoff stage is incremented and the station moves in the state (i + 1, fc), where 
k = 0, . . . , Wi+i — 1, with uniform probability Peq/Wi+i, whereby P^q, i.e., the probability of equivalent failed transmission, 
is defined as P^q = 1 — (1 — Pe)(l — Pcoi) = Pcoi + Pe — Pe ■ Pcoi- Probabilities Pcoi and Pe are, respectively, the collision 
and the packet error probabilities related to the s-th station. 

The transition probabilities for the generic s-th station's Markov process in Fig.[T]could be separated as summarized in what 
follows, depending on whether transitions start from standard backoff states or from post-backoff states. 
Backoff state transitions 
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The meaning of the underlined probabilities is as follows. The first equation in ^ states that, at the beginning of each slot 
time, the backoff time is decremented. The second (third) equation accounts for the fact that after a successful transmission, 
the station goes in post-backoff because of an empty (non empty) queue. In both equations, q is used to identify the probabiUty 

^In order to keep the notation concise, we omit the apex s over the probabilities involved in the model. 
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that the queue contains at least a packet waiting for transmission after a time slot, and it will be better defined in Section 
III-DI where the employed traffic model is described. The fourth equation deals with the situation in which the station has 
reached the retry limit and, after a packet transmission, the buffer of the station is empty. In this situation, the station moves 
in the post-backoff stage with an empty queue. The last equation accounts for a scenario similar to the previous one with the 
difference that, after the packet transmission, the queue is not empty. 
Post-backoff state transitions 

Pp,k\P,k+l 
Po,k\P,k+l 

Pp,o\Pfi 

Pp,k\P,0 
Po,k\P,0 
Pl,k\P,0 

The meaning of the underlined probabilities is as follows. The first equation states that the station remains in the post-backoff 
stage because the queue is empty, whereas the second equation accounts for a transition in the zero-th backoff stage because a 
new packet arrives at the end of a backoff slot. The third equation models the situation in which there are no packets waiting 
for transmission, and the station remains in the state (P, 0) (idle state). 

The fourth equation deals with the situation in which the station is in the idle state {P, 0), and, at the end of a backoff slot, 
a new packet arrives in the queue. In this scenario, the packet is successfully transmitted, the queue is empty, and the station 
moves in another post-backoff stage. The term Pi identifies the probability that the channel is idle, and it is defined as follows 
with respect to the s-th tagged station: 

N 

pt^= n (1--.) 

The fifth equation accounts for a scenario similar to the previous one, except that the station queue is not empty after the 
immediate transmission of a packet or a situation of busy channel. The last equation models the scenario in which the station 
goes from the idle state {P, 0) to the first backoff stage because of a failure of the immediate transmission of the packet arrived 
in the head of the queue. 

B. Throughput Evaluation 

Next line of pursuit consists in finding the probability Tg that the s-th station starts a transmission in a randomly chosen time 
slot. Due to the lengthy algebra involved in the derivations needed for solving the bidimensional Markov chain, the relation 
that defines Tg has been derived in an additional document available at [40], whereas for conciseness we show the final formula 
in (O (shown at the bottom of this page), along with the other key probabilities needed in this paper. Given in ((Sj, we can 
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evaluate the aggregate throughput S as follows: 

N N 



1 1 
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whereby Tqi, is the expected time per slot, PL'^^^ is the packet size of the s-th station, and pj'*^ is the probability of successful 
packet transmission of the s-th station: 
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The evaluation of the aggregate throughput in (HJl requires the knowledge of the expected time per slot, Tav Its evaluation is 
the focus of the next section. 

C. Evaluation of the Expected Time per Slot 

The expected time per slot, Tav, can be evaluated by weighting the times spent by a station in a particular state with the 
probability of being in that state. First of all, we observe that there are four different kinds of time slots, with four different 
average durations: 

> the idle slot, in which no station is transmitting over the channel, with average duration T/; 

« the collision slot, in which more than one station is attempting to gain access to the channel, with average duration Tc', 
« the slot due to erroneous transmissions because of imperfect channel conditions, with average duration Te; 

> the successful transmission slot, with average duration T5. 

The expected time per slot, Tav, can be evaluated by adding the four expected slot durations: 

Tav =Ti + Tc + Ts + Te. (6) 

We will now evaluate T/, Tc, Ts, and Te- 

Upon identifying with a an idle slot duration, and defining with Ptr the probability that the channel is busy in a slot 
because at least one station is transmitting: 

JV 

Ptr^1-Y[{1~Ts) (7) 

a = l 

the average idle slot duration can be evaluated as follows: 

= (1 - Ptr) ■ o (8) 
The average slot duration of a successful transmission, Tg, can be found upon averaging the probability P'f^ that only the 
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s-th tagged station is successfully transmitting over the channel, times the duration of a successful transmission from the 
s-th station: 

JV 

Ts = Y^ Pi"' (1 - Pi')) ■ T^"' (9) 

s = l 

Notice that the term (1 — P'f^) accounts for the probability of packet transmission without channel induced errors. 
Analogously, the average duration of the slot due to erroneous transmissions can be evaluated as follows: 

JV 

Te^Y^ Pi"' ■ Pi"' ■ T^"' (10) 

s = l 

Let us focus on the evaluation of the expected collision slot, Tc- There are different values of the collision probability 
P^^ , depending on the class of the tagged station identified by d. We assume that in a collision of duration T^'' (class-d 



collisions), only the stations belonging to the same class, or to higher classes (i.e., stations whose channel occupancy is lower 
than the one of stations belonging to the tagged station indexed by d) might be involved. 

In order to identify the collision probability Pq \ let us first define the following three transmission probabilities (P^j^^ , P!^l^\ 
P^R^) under the hypothesis that the tagged station belongs to the class d. Probability Pyjj'*^ represents the probability that at 
least another station belonging to a lower class transmits, and it can be evaluated as 

i— 1 sGn(z) 

Probability Pl^lf' is the probability that at least one station belonging to a higher class transmits, and it can be evaluated as 

^tT=i- n n (1--=) (12) 

i — d-\-l s£n{i) 

Probability P^j^^ represents the probability that at least a station in the same class d transmits: 

p?k'^=i~ n (1-^=) (13) 

Therefore, the collision probability for a generic class d takes into account only collisions between at least one station of class 
d and at least one station within the same class (internal collisions) or belonging to higher class (external collisions). Hence, 
the total collision probability can be evaluated as: 



whereby 



Pc^'^ = (l_-P^i''^)-(l~P^i'^V (15) 

sE7i{d) j^n(d) .j^ s 

represents the internal collisions between at least two stations within the same class d, while the remaining are silent, and 

concerns to the external collisions with at least one station of class higher than d. 
Finally, the expected duration of a collision slot is: 

Tc = Yl ' ■ (17) 

Constant time durations T^'\ T^f^ and T^"*' are defined in a manner similar to [22] with the slight difference that the first 
two durations are associated to a generic station s, while the latter is associated to each duration class, which depends on the 
combination of both payload length and data rate of the station of class d. 

D. Traffic Model 

The employed traffic model for each station assumes a Poisson distributed packet arrival process, whereby the inter-arrival 
times among packets are exponentially distributed with mean 1/Af. In order to greatly simplify the analysis, we consider small 
queue, as proposed in [12], even though the proposed analysis may be easily extended to queues with any length. The traffic of 
each station is accounted for within the Markov model by employing a probabilitj^ q, that accounts for the scenario whereby 

''a superscript (t) is used for discerning tlie probability q among the stations. 



at least one packet is available in the queue at the end of a slot. In our setting, each station is characterized by its own traffic, 
and the probability (7^*^ of the t-th station can be evaluated by averaging over the four types of time slots, namely idle, success, 
colUsion, and channel error time slot. Upon noticing that, with the underlined packet model, the probability of having at least 
one packet arrival during time T is equal to 1 — e^^*'^, g^^*^ can be evaluated as: 

qW = (i_Pj,^).(l_e-A*-)+ 

+ EtiP^'^ (l„pi^)).(l-e-A.-T<-")+ 



(18) 



whereby the probabiUties Ptr, Ps'^\ and P^"^ are, respectively, as defined in (|7]i, dU, and (fT4l i. whereas pj'*'' is the packet 



c 

error rate of the s-th station. 

III. Evaluating the Network Loading Conditions 

In a previous paper [14], we proved that the behaviour of the aggregate throughput in a network of N homogeneous 
contending stations is a linear function of the packet arrival rate A with a slope depending on both the number of contending 
stations and the average payload length. We also derived the interval of validity of the proposed model by showing the presence 
of a critical A, above which all the stations begin operating in saturated traffic conditions. 

This kind of behaviour, with appropriate generalizations, is also observed when multirate and variable loaded stations are 
present in the network. We have to identify a set of conditions for a network to be considered as loaded. We notice in passing 
that this framework is not generally considered in the literature, since most papers assume saturated traffic conditions. A key 
observation from the analysis developed in this section is that in an unloaded network there is no need to guarantee fairness, 
since each contending stations can transmit its packets at its own pace, regardless of its minimum CW, as well as other network 
parameters. 

Under the traffic model described in section III-DI we define unloaded a network in which each contending station has a 
packet rate Af less than or equal to its packet service rate fig'': 

At <4*\ VtG {!,..., iV} (19) 

The reason is simple: this condition ensures that the average packet inter-arrival time is greater than or equal to the average 
service time of the t-th station. In such a scenario, the probabilities of collisions among stations are very low, and each 
contending station is able, on the average, to gain the access to the channel as soon as a new packet arrives in its queue. Notice 
that p.g'' only depends on the packet rates A^ of the other iV — 1 stations other than the tagged one. 

The evaluation of the packet service rate /i^^ in a multirate and heterogeneous network is quite difficult [12] since packet 
arrivals may occur during the stage of post-backoff, as well as during the usual backoff stages accomplished by each station 
before gaining the channel for transmission. Since we are interested in a threshold which differentiates the unsaturated from the 
saturated loading conditions of the stations, we can employ an upper bound defined by the saturation service rate, identified 
as in place of the actual service rate The advantage relies on the observation that such a bound is always evaluated 
considering a post-backoff stage. Indeed, after a packet transmission, a new packet is always available in the queue assuming 

(s) 

saturated traffic; therefore, the service time starts from a post-backoff phase whereby the contention window is Wq . We 

''By homogeneous we simply mean that the network is characteiized by Af stations transmitting with the same bit rate (no multirate hypothesis) and the 
same load. 



notice that the saturation service time always includes the post-backoff stage, thus its duration is longer than the actual service 
time evaluated without considering the post-backoff time. 

Hence, in the remaining part of this section, we evaluate the saturation service rate /i^*^ = 1/Tselv, i.e., the rate at which 
packets are taken from the queue of the t-th station under saturated conditions. 

Upon considering the tagged station identified by the index t e {1, • • ■ , N}, the saturation service time TsUv can be defined 
as follows [28]: 



\ 3=0 



DROP 

The first term in the summation represents the average time that a station spends through the backoff stages from to r before 
transmitting a packet, i.e., the so called MAC access time. We notice that for the i-th stage, i collisions of average duration 
Tc, as well as i backoff stages from to z (each of them with an average number Wj*' of slot of duration TJ^'*^) occurred, 
after which the packet is successfully transmitted with duration t]*''. The second term of the summation takes into account 
the average duration of a packet drop that occurs after (r + 1) collisions and backoff stages. The whole summation is scaled 
by a normalization factor that takes into account the probability set over which the service time is evaluated. 
The average number of slots for the i backoff stages, is defined as 

Wf^ = (2™nO-.™) . w^*'> - l)/2. 

Each slot has average duration T)|*\ which is substantially evaluated as Tav in (ID except that the tagged station it) is not 
considered because it is either idle, or in a backoff state. 

Let us discuss two sample scenarios in order to derive a variety of observations that are at the very basis of the fairness 
problem developed in the next section. The network parameters used in the investigated IEEE802.11b MAC layer are reported 
in Table |l][l]. The first investigated scenario considers a network of 3 contending stations. Two stations, namely S2 and S3, 
transmit packets with constant rates A = 100 pkt/s and A = 500 pkt/s, respectively. The bit rate of the two stations S2 and 
S3 is 11 Mbps. The third station, SI, has a bit rate equal to 1 Mbps and a packet rate Ai that is varied in the range [0, 2000] 
pkt/s in order to investigate its effects on the network load. The behaviour of the three ratios Xt/f^"P is shown in Fig.|2l Some 

(2) (3) 

observations are in order First of all, notice that as far as S 1 increases its packet rate, the curves erv and A3 ■ Tserv tend 

to increase because of the increasing values of the service times Tserv and Tserv experienced by S2 and S3. Notice that, as 
Ai increases, the slowest station SI tends to transmit more often. The station S3 goes in loaded condition when Ai « 20pkt/s, 
while S2 can be considered loaded for Ai ~ 200 pkt/s. When a station becomes loaded, the incoming packets tend to be stored 
in the station queue waiting for transmission since the service rate, i.e., the number of packets that on average are serviced by 
the MAC, is below the rate by which the packets arrive in the station queue. 

The per-station throughput achieved by the three stations in the investigated scenario is shown in Fig. |3] The throughput 
gained by the two fastest stations, S2 and S3, tends to decrease because of the anomaly problem in the multirate scenario: 
the slowest station tends to occupy the channel longer and longer as far as its packet rate Ai increases. In the same figure, 
we show two tick curves. The horizontal line L2 corresponds to the saturation throughput of SI, while the straight line LI is 
the tangent to the throughput curve passing through the origin. For very small values of Ai, the throughput of the station SI 
grows linearly with Ai. Packets are mainly transmitted as soon as they arrive at the MAC layer, and the station throughput is 



approximately equal to Ai ■ PL^^\ However, when the station approaches the transition point A^; = 89.19 pkt/s derived with the 
proposed framework (this is the value of Ai corresponding to the relation Ai//i^^^ — 1), the throughput curve tends to reach 
the asymptote L2, which corresponds to the saturation throughput of S 1 . Notice that the curve L2 approximately corresponds 
to 0.73 Mbps, which is A| • PL^^^ = 89.19 • 1028 • 8 bps. 

In the second scenario, the stations SI and S2 are interested by a constant packet rate noticed in the label of Fig. |4] These 
two stations do not get loaded by the increasing packet rate of the station S3 since the curves \i/ ji'g^ and X2/n^P are strictly 
less than one. As a consequence, the per-station throughput of both S 1 and S2 is approximately constant across the range of 
values of the packet rate A3 as noticed in Fig. |5] On the other hand, the throughput achieved by the station S3 tends to saturate 
as soon as A3 reaches the value 533 pkt/s noticed in Fig. |4l 

In the light of the previous two sample scenario, let us summarize the main ingredients of the results proposed in this 
section. As observed in the two previous sample scenarios, this method allows to identify whether the network is loaded by 
estabUshing the thresholds of each contending station in the network. This issue has been overlooked in the literature, where 
the saturation assumption is widely adopted. Moreover, this issue is at the very basis of any throughput optimization strategy 
since an unloaded network does not need to be optimized. 

We say that the network is loaded when every station has a traffic above its proper threshold. On the other hand, it should 
be noticed that a network can be unloaded even if a subset of the stations is loaded. This was the case of the second scenario 
described above, where, despite the fact that the station S3 was interested by an increasing traffic load A3, the stations SI and 
S2 did not experience any performance loss (see Fig. |5]l because their traffics were below the respective thresholds. 

IV. The Proportional Fairness Throughput Allocation Algorithm 

This section presents the novel throughput allocation criterion, which aims at improving fairness among the N contending 
stations. In order to face the fairness problem in the most general scenario, i.e., multirate DCF and general station loading 
conditions, we propose a novel Proportional Fairness Criterion (PFC) by starting from the PCF defined by Kelly in [41], and 
employed in [25] in connection to proportional fairness throughput allocation in multirate and saturated DCF operations. 

It is known that one of the main drawbacks of the basic DCF operating in a multirate scenario relies on the fact that it 
behaves in such a way as to guarantee equal long-term channel access probability to the various contending stations [25], 
[30]. In order to solve this problem, various optimization algorithms have been proposed in the literature (see, for instance, 
[25]-[37]). These contributions allowed to highlight the behaviour of the DCF as well as various drawbacks when operating in a 
multirate scenario. For instance, it is known that the aggregate throughput of multirate IEEE802.1 1-like networks is optimized 
when only the high rate stations transmit, while the low rate stations are kept silent. Of course, this result is not desirable from 
a fairness point of view, even though it mitigates the DCF performance anomaly noticed in [2]. 

To the best of our knowledge, the solutions proposed so far in the literature refer to homogeneous networks in that all 
the contending stations operate with the same traffic, as exemplified by the station packet rate A^, and channel conditions. 
Furthermore, almost all the works focus on saturated traffic conditions. As already mentioned before, in a practical setting the 
contending stations have their own traffic and are affected by different channel conditions. The key observation here is that a 
fair throughput allocation should account for the station packet rate, as well as for the specific channel conditions experienced 
by each contending station. 

Let us briefly mention the rationales at the very basis of the PFC. A proportional fairness optimization criterion allocates to 
each station a throughput proportional to the station transmission rate. Resorting to the notation proposed in [41], a throughput 



allocation vector x — {xg', s — 1, - ■ ■ , N} is proportional fair if the following condition holds: 



N ^ _ 

<0 (21) 



-I -^s 
s— 1 



for any other feasible throughput allocation vector y*. The PFC maximization problem, which satisfies S2T[ . can be formalized 
as follows: 

N 

max ^log(a;3) (22) 

s=l 

over Xs e[0,Xs^m], s = l,...,N 

whereby Xg^m is the maximum throughput of the s-th station. Due to the strict concavity of the logarithmic function and because 
of the compactness of the feasible region Xg £ [0,Xs.m], s ~ 1, . . . ,N, there exists a unique solution to the optimization 
problem (l22T i. This implies that a local maximum is also global. 

Given these preliminaries, in the proposed model, the traffic of each station is characterized by the packet arrival rate As, 
which depends mainly on the application layer Let Xmax be the maximum packet rate among Ai, . . . , Xn- 

Consider the following modified optimization problem: 

u^uiSu--- ,SN)^j:tijt::-^ogiSg) ^^^^ 

subject to Ss e [0, Ss^rn], S^1,...,N 

(s) 

whereby Sg is the throughput of the s-th station, and Ss.m is its maximum value, which equals the station bit rate -R^ . In 
our scenario, the individual throughputs, Sg, are interlaced because of the interdependence of the probabilities involved in the 
transmission probabilities r^, Vs — 1, . . . ,N. For this reason, we reformulate the maximization problem in order to find the 
N optimal values of Ts for which the cost function U in (l23T l gets maximized. The optimal values r* are then used to set the 
network parameters of each station. 

Due to the compactness of the feasible region Sg G [0, S'^.-m], Vs, the maximum of U{Si, ■ ■ ■ , Sn) can be found among the 
solutions of \7U = (^^, ■ ■ ■ , = 0- After some algebra (the derivations are reported in the Appendix), the solutions can 

be written as: 

' E T^ = ^f^' V,^l,...,iV (24) 



whereby C — X^iLi ^^'^ T^av is a function of ri, • ■ • ,tn as noticed in (O. 

Due to the presence of Tav, a closed form of the maximum of U{Si, ■ ■ ■ ,Sn) cannot be found. Notice that it is quite 
difficult to derive the contribution of the partial derivative of Tav on Tj, especially when ^ 1, because of the huge number 
of network parameters belonging to different stations. The definition of Tav in I© is composed by four different terms, which 
include the whole set of t^. Vs. In order to overcome this problem, we first numerically obtain the optimal values t*, Vs from 
(l24l i. then, we choose the value of the minimum contention window size, VFq"'', by equating the optimizing t* to ^ for any 
s. 

The optimization criterion summarized in ( l23T l will be denoted as Load Proportional Fairness (LPF) criterion in the following. 

Let us derive some observations on the proposed throughput allocation algorithm by contrasting it to the classical PF 
algorithm. Consider two contending stations with packet rates Ai = 50 pkt/s and A2 = 100 pkt/s, respectively. Employing 
the classical PF method, a throughput allocation is proportionally fair if a reduction of x% of the throughput allocated to one 
station is counterbalanced by an increase of more than x% of the throughputs allocated to the other contending stations. 



In our setup, the ratio A1/A2 can be interpreted as the frequency by which the first station tries to get access to the channel 
relative to the other station. Therefore, a throughput allocation is proportionally fair if, for instance, a reduction of 20% of the 
throughput allocated to the first station, which has a relative frequency of 1/2, is counterbalanced by an increase of more than 
40% of the throughput allocated to the second station. In a scenario with multiple contending stations, the relative frequency 
is evaluated with respect to the station with the highest packet rate in the network, that gets unitary relative frequency. 

Based on extensive analysis, we found that the optimization problem (1231 ) sometimes yields throughput allocations that cannot 
be actually managed by the stations. As a reference example, assume that, due to the specific channel conditions experienced, 
the first station has a bit rate equal to 1 Mbps and needs to transmits 200 pkt/s. Given a packet size of 1024 bytes, that is 
8192 bits, the first station would need to transmit 8192 x 200 bps « 1.64Mbps far above the maximum bit rate decided at 
the physical layer In this scenario, such a station could not send over the channel a throughput greater than 1Mbps. The 
same applies to the other contending stations in the network experiencing similar conditions. In order to face this issue, we 
considered the following optimization problem 



and A^Qj. = maxs A*. The allocation problem in (IZST i, solved as for the LPF in (|23] |. guarantees a throughput allocation which 
is proportional to the frequency of channel access of each station relative to their actual ability in managing such traffic. The 
idea relies on the observation that it is not optimal to allocate a throughput that the station cannot actually manage. 

The optimization criterion summarized in dZSl l will be denoted as Modified Load Proportional Fairness (MLFF) criterion 
in what follows. As in the case of the LPF criterion, the MLPF optimization problem in dZSt is solved by first numerically 
obtaining the optimal values t*.Vs, and then choosing the value of the minimum contention window sizes Wq'''' by equating 
the optimizing t* to (O independently for any s. 



This section presents simulation results obtained for a variety of network scenarios optimized with the fairness criteria 
proposed in the previous section. 

We have developed a C++ simulator modelling both the DCF protocol details in 802.11b and the backoff procedures of a 
specific number of independent transmitting stations. The simulator considers an Infrastructure BSS (Basic Service Set) with 
an Access Point (AP) and a certain number of fix stations which communicate only with the AP. Traffic is generated following 
the exponential distribution for the packet interarrival times. Moreover, the MAC layer is managed by a state machine which 
follows the main directives specified in the standard [1], namely waiting times (DIFS, SIFS, EIFS), post-backoff, backoff, basic 
and RTS/CTS access modes. The typical MAC layer parameters for IEEE802.11b reported in Table |l][l] have been used for 
performance validation. 

The first investigated scenario, namely scenario A, considers a network with 3 contending stations. Two stations transmit 
packets with rate A = 500 pkt/s at 11 Mbps. The payload size, assumed to be common to all the stations, is PL = 1028 
bytes. The third station has a bit rate equal to 1Mbps and a packet rate A — 1000 pkt/s. 



max Y.s=i ■ ^og{Ss) 



(25) 



over 5*, e [0,Ss,m], s = l,...,N 



where Vs = 1, . . . , 



N, 




V. Simulation Results 



From ( fT9b . it is straightforward to notice that the scenario A refers to a loaded network. Fig. |6] shows the service rates /i^^ 
of the three contending stations as a function of the packet rate Ai of the station transmitting at 1 Mbps. The operating point 
of the considered scenario A is highlighted in Fig. |6] Notice that, since the service rates /i^*'' of the three stations are below 
the respective packet rates Af, the network is loaded. Moreover, the service rates of the three stations tend to the same values 
because of the rate anomaly problem: the station transmitting at 1 Mbps reduces the service rates of the other stations. 

The simulated normalized throughput achieved by each station in this scenario is depicted in the left subplot of Fig.|2]for the 
following four setups. The three bars labelled 1-DCF represent the normalized throughput achieved by the three stations with 
a classical DCF The second set of bars, labelled 2-PF, identifies the simulated normalized throughput achieved by the DCF 
optimized with the PF criterion [25], [41], whereby the actual packet rates of the stations are not considered. The third set 
of bars, labelled 3-LPF, represents the normalized throughput achieved by the three stations when the allocation problem (l23T l 
is employed. Finally, the last set of bars, labelled 4-MLPF, represents the simulated normalized throughput achieved by the 
contending stations when the CW sizes are optimized with the modified fairness criterion in ( l25T l. Notice that the throughput 
allocations guaranteed by LPF and MLPF improve over the classical DCF. When the station packet rate is considered in 
the optimization framework, a higher throughput is allocated to the first station presenting the maximum value of A among 
the considered stations. However, the highest aggregate throughput is achieved when the allocation is accomplished with the 
optimization framework 4-MLPF. The reason for this behaviour is related to the fact that the first station requires a traffic 
equal to 8.22 Mbps — 10^ pkt/s • 1028 bytes/pkt • 8 bits/pkt, which is far above the maximum traffic (1 Mbps) that the station 
would be able to deal with in the best scenario. In this respect, the MLPF criterion results in better throughput allocations 
since it accounts for the real traffic that the contending station would be able to deal with in the specific scenario at hand. 

Similar considerations can be drawn from the results shown in the right subplot of Fig. [7] (related to scenario B), whereby 
in the simulated scenario the two fastest stations are also characterized by a packet rate greater than the one of the slowest 
station. Notice that the optimization framework 3-LPF is able to guarantee improved aggregate throughput with respect to both 
the non-optimized DCF and the classical PF algorithms. The operating point of the scenario B is highlighted in Fig. |6l based 
on the considerations above, this is a loaded network as well. 

The aggregate throughputs achieved in the two investigated scenarios are as follows: 



Scenarios 


1-DCF 


2-PF 


3-LPF 


4-MLPF 


A 


Jain's Index 
S [Mbps] 


0.451 
1.85 


0.872 
3.60 


0.724 
3.10 


0.848 
5.07 


B 


Jain's Index 
S [Mbps] 


0.474 
1.99 


0.881 
3.63 


0.976 
4.62 


0.847 
5.07 



whereby we also show the fairness Jain's index [42] evaluated on the normalized throughputs noted in the subplots of Fig.|71 It 
is worth noticing that the proposed MLPF throughput allocation criterion is able to guarantee improved aggregate throughput 
relative to both the classical DCF and the PF algorithm, with fairness levels on the same order of the ones guaranteed by the 
classical PF algorithm. 

For the sake of investigating the behaviour of the proposed allocation criteria as a function of the packet rate of the slowest 
station, we simulated the throughput allocated to a network composed by three stations, whereby the slowest station, transmitting 



at 1Mbps, presents an increasing packet rate in the range 10—3000 pkt/s. The other two stations transmit packets at the constant 
rate A = 500 pkt/s at 1 1 Mbps. The simulated throughput of the three contending stations is shown in the three subplot of 
Fig. [8] for the unoptimized DCF, as well as for the two criteria LPF and MLPF Some considerations are in order. Let us focus 
on the throughput of the DCF (uppermost subplot in Fig. [8]l. As far as the packet rate of the slowest station increases, the 
throughput allocated to the fastest stations decreases quite fast because of the performance anomaly of the DCF [2]. The three 
stations reach the same throughput when the slowest station presents a packet rate equal to 500 pkt/s, corresponding to the 
one of the other two stations. From A = 500 pkt/s all the way up to 3000 pkt/s, the throughput of the three stations do not 
change anymore, since all the stations have a throughput imposed by the slowest station in the network. Let us focus on the 
results shown in the other two subplots of Fig. [8] labelled LPF and MLPF, respectively. A quick comparison among these three 
subplots in Fig. [8] reveals that the MLPF allocation criterion guarantees improved aggregate throughput for a wide range of 
packet rates of the slowest station, greatly mitigating the rate anomaly problem of the classical DCF operating in a multirate 
setting. In terms of aggregate throughput, the best solution is achieved with the MLPF criterion, which avoids that the slowest 
station receives a throughput allocation that would not be able to employ due to its reduced bit rate (1 Mbps). 

VI. Conclusions 

Focusing on multirate IEEE 802.11 Wireless LAN employing the mandatory Distributed Coordination Function (DCF) 
option, this paper established the conditions under which a network constituted by a certain number of stations transmitting 
with their own bit rates and packet rates can be considered loaded. It then proposed a modified proportional fairness criterion 
suitable for mitigating the rate anomaly problem of multirate loaded IEEE 802.11 Wireless LANs. 

Simulation results were presented for some sample scenarios showing that the proposed throughput allocation was able to 
greatly increase the aggregate throughput of the DCF, while ensuring fairness levels among the stations of the same order of 
the ones available with the classical proportional fairness criterion. 
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Appendix 

The objective of this section is to derive the maximum of the function U{Si^ • • ■ , Sm) in ( |23] ) as a function of Ts, Vs 
1,...,7V. 

Upon substituting Q and I© in (|4|l, and deriving with respect to Tj, we obtain: 
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(26) 



whereby the last relation stems from the independence of Pe and PL^*) on tj. Exchanging the derivative with the summation 
yields; 

Aj 1 1 '^k C dTav w ■ 1 Ar /'^-7^ 



whereby C — x^^ ■ equating (l27l i to zero, (|24] | easily follows. 
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Fig. 1. Markov chain for the contention model of the generic s-th station in general traffic conditions, based on the 2- way handshaking technique, considering 
the effects of channel induced errors, unloaded traffic conditions, and post-backoff. 
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Fig. 2. Behaviour of the ratio At//jg' = At ■ T^^Xv in 1191 in a network of three contending stations (labelled SI, S2 and S3) as a function of the packet 
rate Ai of the slowest station SI transmitting at 1 Mbps. The other two stations transmit at 11 Mbps with constant packet rates, respectively equal to lOOpkt/s 
and 500pkt/s. The packet size PL is equal to 1028 bytes for the three contending stations. 




Fig. 3. Behaviour of the per-station throughput in a network of three contending stations (labelled SI, S2 and S3) as a function of the packet rate Ai of the 
slowest station SI transmitting at 1 Mbps. The other two stations transmit at 11 Mbps with constant packet rates, respectively equal to lOOpkt/s and 500pkt/s. 
The packet size PL is equal to 1028 bytes for the three contending stations. 
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Fig. 4. Behaviour of the ratio At//ig' = \t - TsUy in il9\ in a network of three contending stations (labelled SI, S2 and S3) as a function of the packet rate 
A3 of the fastest station S3 transmitting at 11 Mbps. The station SI transmits with constant packet rate lOpkt/s at 1 Mbps, whereas the station S2 transmits 
with constant packet rate lOOpkt/s at 11 Mbps. The packet size PL is equal to 1028 bytes for the three contending stations. 




Fig. 5. Behaviour of the per-station throughput in a network of three contending stations (labelled SI, S2 and S3) as a function of the packet rate A3 of the 
fastest station S3 transmitting at 11 Mbps. The other two stations transmit with constant packet rates, respectively equal to lOpkt/s and lOOpkt/s. The packet 
size PL is equal to 1028 bytes for the three contending stations. 



TABLE I 



Typical network parameters 



MAC header 


28 bytes 


Propag. delay Tp 


1 /is 


PLCP Preamble 


144 bit 


PLCP Header 


48 bit 


PHY header 


24 bytes 


Slot time 


20 /iS 


PLCP rate 


1Mbps 


Wo 


32 


No. back-off stages, m 


5 




1024 


Payload size 


1028 bytes 


SIFS 


10 us 


ACK 


14 bytes 


DIFS 


50 fis 


ACK timeout 


364/iS 


EIFS 


364 fj,s 




100 200 250 300 400 500 600 700 800 900 1000 

[pkt/s] 



Fig. 6. Behaviour of the critical packet rates (/x^' = l/r|g^„) in a network of three contending stations as a function of the packet rate Ai of the slowest 
station transmitting at 1 Mbps. The other two stations transmit with a constant packet rate equal to 500pkt/s at 1 1 Mbps. Notice that curves related to 
the stations at 1 1 Mbps are superimposed since they both employ the same network parameters. 
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Fig. 7. Simulated normalized thi'oughput achieved by three contending stations upon employing 1) a classical DCF; 2) DCF with PF 
optimized as noted in )23L and 4) DCF optimized with the MLPF criterion. Left and light plots refer to scenarios A and B, respectively. 
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Fig. 8. Simulated normalized throughput achieved by three contending stations as a function of the packet rate of the slowest station in DCF, LPF and MLPF 
modes. 



